o 

(N 

> 
O 

(N 



X 

6 

o3 



> 
O 

in 
o 



X 



32nd International Cosmic Ray Conference, Beijing 201 1 
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Abstract: The bulk of air showers initiated by very high energy cosmic rays exhibits a longitudinal development in depth 
with a single well-defined shower maximum. However, a small fraction of showers has a profile that differs considerably 
from this average behaviour. In extreme cases, such anomalous longitudinal profiles can even have two distinct shower 
maxima. We discuss the properties of the primary interactions that lead to such profiles. Simulations are used to estimate 
the rate of anomalous profiles in dependence of primary energy, mass, and characteristic features of hadronic multiparticle 
production at very high energies. 
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1 Introduction 

Cosmic ray detectors like HiRes [1], the Pierre Auger Ob- 
servatory [2], or Telescope Array [3] measure the longi- 
tudinal profiles of air showers initiated by ultra-high en- 
ergy cosmic ray particles. The data of these experiments 
offers the possibility to study hadronic interactions beyond 
the reach of man-made accelerators. In practice, however, 
the interpretation of the data is difficult because of the un- 
known mass composition of the cosmic particle beam. 

In this paper, we present a Monte Carlo study of longi- 
tudinal profiles using the CONEX air shower simulation 
program [4] to study if 'anomalous' longitudinal shower 
profiles could provide a unique signature for properties of 
hadronic interactions. Measurements [5], simulations [6] 
and theoretical considerations (e.g. [7]) suggest that the 
longitudinal development of the electromagnetic compo- 
nent of ultra-high energy air showers is to a good approxi- 
mation universal, i.e. that the shape of the profiles is very 
similar for the bulk of showers irrespective of the primary 
energy, particle type or hadronic interaction model. In rare 
cases, however, simulated profiles can significantly deviate 
from this universal shape and anomalous shower profiles 
can be found that in extreme cases show two completely 
separated maxima. 

Examples of such anomalous shower profiles are shown 
in Fig. 1. The profile in the left panel is generated by a 
primary helium nucleus with an energy of 10 17 eV. The 
first interaction occurs at 19 g cm~ 2 and results in a sub- 
shower with a maximum at around 700 g cm~ 2 . A specta- 
tor nucleon with a quarter of the primary energy penetrates 
deeply into the atmosphere before it interacts at 583 g cm -2 
and creates a second sub-shower that reaches its maximum 



beyond lOOOg cm -2 . Similar cases are shown in Fig. 1(b) 
and (c), but now for primary protons. Here, instead of a 
spectator nucleon, the second shower is created by the lead- 
ing particle from the first interaction. 

These examples suggest, that the anomalous shower pro- 
files originate from leading or spectator particles that pen- 
etrate deeply into the atmosphere before interacting. The 
probability for propagating a slant depth distance greater 
than AX without interacting is given by 



P(AX) = e 



= P A 



(1) 



where A denotes the hadronic interaction length in air. As- 
suming that a propagation distance of at least 300 g cm~ 2 
is needed to experimentally distinguish these showers from 
the universal shower profile, one would expect P ~ 3-10 -3 
and P ~ 8 ■ 10~ 4 for 10 17 eV and 10 19 eV primaries using 
a hadronic interaction length as predicted by the SIBYLL 
interaction model. The actual rate of occurrence will how- 
ever be smaller, because for an experimental observation of 
anomalous shower profiles, the inelasticity of the interac- 
tion must be in a suitable range such that both sub-showers 
carry a fraction of the primary energy that is large enough 
for a detection. 

2 Analysis of Air Shower Simulations 

To estimate the frequency of detectable anomalous air 
showers we generated more than 10 4 showers for each 
combination of three different hadronic interaction models 
(Sibyll2.1 [8], QgsjetII [9] and Epos 1.99 [10]), seven 
primary energies between E = 10 17 and 10 20 eV and three 
primary masses (proton, helium, and iron). 
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(a) Helium E = 10 17 eV (EPOS) 




X in g/cm 2 

(b) Proton E = 10 19 5 eV (SIBYLL) 




X in g/cm 2 

(c) Proton E = 10 20 eV (SIBYLL) 



Figure 1 : Examples of anomalous shower profiles simulated with CONEX. 



Anomalous profiles are identified by comparing the \ 2 - 
value obtained from a fit of the simulated profile with a 
Gaisser-Hillas function [11], /gh, with the x 2 that results 
from the sum of two Gaisser-Hillas functions: 

2 _ ^ (iVi-/GHl(*i)-/GH2(*«)) 2 n . 

Xdouble — 2-^1 t7 ' ' 

1 = 

with 

/gHI = /GH (-^rnax,l! ^max,lj Ai, Xq) (3) 
/GH2 = /GH (^V ma x,2i AX max , A 2 , X 2 ) , (4) 

where Ni denotes the simulated shower size at slant depth 
Xi and n is the number of generated data points. The nor- 
malisation of the functions is given by -ZV max i and X max l 
is the depth of maximum of the first sub-shower. Follow- 
ing [12], the shape parameters Xi and Xq are constrained 
to average values corresponding to the 'normal' universal 
profiles to assure a better stability of the fit. Finally, AX max 
denotes the difference between the shower maxima of the 
first and second sub-shower. 

For an anomalous profile that originates from a deeply pen- 
etrating leading particle, we expect that the fit with two 
Gaisser-Hillas functions will lead to parameters that are re- 
lated to the properties of the two sub-showers. The inelas- 
ticity k of the first interaction can be estimated from the 
normalisations of the two Gaisser-Hillas functions, 

~ 1 -^max,2 

^ TV TN ' (5) 

J >max.l T J»max.2 

and the fitted distance between the two shower maxima is 
related to the slant depth distance between the first and sec- 
ond interaction 

AX » AJW (6) 

This relations assumes that the shower development from 
the interaction point is a mere translation, that is the same 
for both showers. This should be a good approximation, 
since the average slant depth distance between the interac- 
tion point and the shower maximum depends only logarith- 
mically on energy. 



In order to obtain a meaningful \ 2 m Eq. (2), a choice 
must be made for V{. CONEX generates 'exact' shower 
sizes without statistical fluctuations. A Poissonian ansatz, 
Vi = Ni would lead to an increasing sensitivity of the x 2 s 
to deviations from a universal profile with energy. For this 
study however, we prefer to have a constant sensitivity 
at each energy and therefore use Vi = kNi/E, where k 
is chosen such that VE^/E N i = °- 01 - with this 
re-normalisation, equal \ 2 probabilities for all energies 
are assured that enable us to study the 'intrinsic' evolution 
of the fraction of anomalous profiles with energy. Fur- 
thermore, profile points with \/V~i/Ni > 0.3 are excluded 
from the fit to avoid the influence of the long tail of 
deeply penetrating muonic component of the shower size 
profile that can not be well described by a four parameter 
Gaisser-Hillas function; Fluorescence detectors can not 
observe this tail above the night-sky background. 

For the analysis of the generated showers, each profile 
point is fluctuated by y/Vl using Gaussian random numbers 
and the A\ 2 between single- and double-Gaisser-Hillas fit 
is calculated. Then we define a shower as 'anomalous' if 

• two Gaisser-Hillas functions significantly improve 
the goodness of fit (A^ 2 > 25), 

• the shower maxima are clearly separated (AX max > 
300 gem" 2 ), 

• both fitted X m . dx were found within the profile range, 

• both fitted sub-showers have more than 20% of the 
primary energy. 

3 Results 

Using the criteria from the last section to select and define 
anomalous shower profiles, their fraction of the total simu- 
lated sample can be determined. As can be seen in Fig. 2, 
no anomalous showers are predicted for iron nuclei, which 
is not surprising because spectator nucleons can carry less 
than 2% (1/56) of the primary energy. Heavier fragment 
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(a) SIBYLL (b) QGSJETlI (c) EPOS 

Figure 2: Fraction of anomalous shower profiles for different primaries and interaction models. 



nuclei could have a larger energy, but due to their large 
cross section they have a small mean free path length in the 
atmosphere. 

Anomalous shower profiles of helium and proton primaries 
occur with a very similar frequency and the fractions 
are of the same order of magnitude as estimated with 
the simple analytic model, Eq. 1. The fraction decreases 
with energy as it would be expected due to the rise of 
the cross sections with energy in the simulations. In case 
of SIBYLL and EPOS, the fraction rises slightly again at 
the highest energies. The reason for this behaviour could 
be the Landau-Pomeranchuk-Migdal (LPM) effect [13] 
on photons produced by the decay of a leading ir°. The 
LPM effect significantly increases the mean free path of 
electromagnetic particles above 10 18 — 10 19 eV, and in 
the case of SIBYLL and EPOS a leading ir° can carry more 
than 50% of the primary energy. In the case of QGSJETlI 
no leading ir° are produced. Further studies are needed to 
confirm this hypothesis. 

The correlation between the fitted AX mal and the slant 
depth difference AX = X n — X^ st of the two most rele- 
vant interaction points is shown in Fig. 3. Here X^ st is the 
interaction depth of the primary particle in the atmosphere 
and X n the interaction depth of the nth leading particle 1 , 
where interactions are only considered if they have a mini- 
mal inelasticity of k > 0.15. As can be seen, the distance 
between the shower maxima is a good estimator for AX 
for the majority of the showers. In addition to the events 
populating the diagonal, there is a cluster of events at small 
AX but large AX mal . These events are due to deeply pene- 
trating sub-showers created not by leading particles, so X n 
and thus also AX are not the correct quantities. 

As a further test of the sensitivity of the anomalous profile 
fraction to the hadronic interaction length, the response of 
the event fraction to a decrease of the hadronic cross sec- 
tions in the simulation is studied. For this purpose all cross 
section above 10 15 eV are rescaled logarithmically during 
CONEX simulations by a factor f(E), 



where E is the energy and /ig is the rescaling factor of the 
cross section at 10 19 eV [14]. The resulting fractions for 
diminished cross sections are shown in Fig. 4. As can be 
seen, decreasing the cross-section of SIBYLL to 80% means 
a 50% and 100% higher fraction at 10 18 eV and 10 19 eV 
respectively and a reduction to 60% increases the rate even 
further. 



4 Conclusions 

In this paper we presented a simulation study on anoma- 
lous longitudinal shower profiles. These showers exhibit 
in extreme cases two distinct shower maxima and originate 
from deeply penetrating spectator nucleons or leading par- 
ticles. It was shown that their rate is expected to be largest 
at low energies and for light primary masses. It is however 
worthwhile noting that the absolute rate estimates given in 
Fig. 2 are only indicative and depend on the selection crite- 
ria as well as on the precision of the profiles, Vi. A precise 



1 . For n = 1 this refers to the leading particle of the primary 
interaction, for n — 2 to the leading particle of the interaction of 
the leading particle of the primary interaction, etc.. 




f(E) = 1 + (/ 19 - 1) 



In (E/10 15 eV) 
ln((10 19 eV/10 15 eV)' 



(7) 



100 200 300 400 500 600 700 800 
X n- X first g/Cm 2 

Figure 3: Correlation between AX and AX max for protons 

(SIBYLL). 
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Figure 4: Anomalous profile fraction for a modified cross 
section in SIBYLL. 

estimate of the rate of anomalous showers detectable by 
current detectors requires to study detailed detector simu- 
lations, which is beyond the scope of this work. 

The experimental detection of such showers would be an 
unambiguous proof of the presence of light nuclei in the 
cosmic particle beam. Moreover, if it is possible to detect 
enough of these events, they could provide a novel tool to 
study the cross section and inelasticity of hadronic interac- 
tions at ultra-high energies. 
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